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ABSTRACT 



A dynamic pressure decaying exponentially with time, 
whose analytic representation in terms of non-dimensional 
quantities is: 

P = P 0 e_T (t > 0) , 

and whose peak value p 0 is larger than the static collapse 
pressure p s , is applied uniformly throughout a spherical 
shell made of rigid-perfectly plastic material. The problem 
is to determine the final deflection of the shell for a 
given yield surface, the shell has no holes and is simply 
supported at the outer edge. 

Three approximate yield surfaces have been used: the 

two-moment limited interaction (HS) , the uncoupled square 
(SS) and the uncoupled diamond (DS) yield surfaces. 

For each yield surface selected, the static collapse 
pressure is first determined both for simply supported and 
clamped edges. In the dynamic problem, the shell is simply 
supported . 

The first part deals with shallow shells. With the 
(HS) and (SS) yield surfaces, solutions for three ranges of 
pressure have been obtained. With the (DS) yield surface, 
solution for one range of pressure has been obtained. 

The second part deals with deeper spherical caps. For 
all three yield surfaces, only solutions for one range of 
pressure have been obtained. 

For shallow shells , graphs giving maximum central 
deflection and energy absorbed as functions of the pressure 
difference p 0 -p s and maximum central deflection as functions 
of energy absorbed are also presented. 
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Marine Engineering 
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NOTATION 



0 




Global and Local Coordinates of the Reference Surface 

Fig . 1 



Oz: Axis of symmetry. 

M: Generic point on the reference surface, which is the 

mid -surface for a homogeneous shell of constant 
thickness . 

Distance from M to the axis of symmetry Oz. 

Distance from M to the tangent plane at the apex 0. 
Unit normal to the reference surface, positive inward, 
Unit tangent to the meridian at M, pointing away from 
the apex 0 . 

Unit tangent to the parallel at M. 

fc 0 = \ x n 

U^: Displacement component in the t^ direction. 

U : Displacement component in the n direction. 



r : 
z : 
n 

<i> : 

v 

* 

n 
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A: 

T : 
o 


A reference length. 

X 

A reference time. 


ii 

> 


Non-dimensional meridional component of the 




displacement . 


u 

n 

« = — : 


Non-dimensional normal component of the displace- 



ment . 



t : 


Time variable. 


_ t 

T “ T : 

o 


Non-dimensional time variable. 


r 

X = A : 


Non-dimensional space variable. 


z 

y = a * 


Non-dimensional space variable. 


dv 

v = : 


Non-dimensional meridional component of the 



velocity . 



9w 

w = : 

3t 


Non-dimensional normal component of the velocity. 


e 6 : 

V 

V 


Non-dimensional strain rate m the tg direction. 

Non-dimensional strain rate in the t^ direction. 

Non-dimensional curvature rate in the tg direction 

->■ f 

Non-dimensional curvature rate in the t^ direction 
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reference 
surface 




Fig. 2a 



Fig . 2b 



Local Stresses 



u e* 

a 

< P 

r . : 
n 4> 

2H : 

5: 



Membrane Stress in the t Q direction. 

y 

Membrane stress in the t, direction 

<f> 

Transverse shear stress. 

Shell thickness. 

Distance from a point of the shell to the reference 
surface, positively inward. 

For thin shell, for which ;r— << 1, where R is the 

R p 

P 

smaller principal radius of curvature of the reference 
surface, the stress and moment resultants per unit 
length are: 



N e - 



f +H 






f +H 




f 


a fl d? 
J-H 9 


/ 


N. = 

$ J 


L„v s 


/ 


II 

CO 



+H 



-H 



G A d ^ 
n4> 



M = 



+H 



-H 



a 0 £d£ 



m a = 
<p 



+H 



J -H 



V d5 



a : Yield stress . 

o 

N : a x 2H 

o o 



M : a H 

o o 
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^int 



ext 



N 



n n = rr— : Non-dimensional stress resultant in the 
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N 



N, 



iN d> 

n, = rr 1 : Non-dimensional stress resultant in the 



* 



s = 



m . = 



N 



N 



M 



M ^> 



m , = 






M 



o 



D . . : 
int 



D 



ext ' 



D 



int 



2ttN A 2 
o 



tg direction. 



t, direction. 

4> 

Non-dimensional transverse shear stress 
resultant . 

Non-dimensional moment resultant of the 
Cg-stresses . 

Non-dimensional moment resultant of the 
a , -stresses . 

<t> 

Rate of internal energy dissipation. 
Rate of external energy input. 

Non-dimensional rate of internal energy 
dissipation . 



D 



ext 



2nN A 2 
o 



T 

o 



Non-dimensional rate of external energy 
input . 



p: Surface density of the shell material. 

P: Applied external pressure, in the normal 

direction. 




s : 



Non-dimensional pressure applied exter- 
nally in the normal direction. 
Non-dimensional applied traction at the 
edge, in the n-direction. 
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n , : 

<P ■ 


Non-dimensional applied traction at the edge, 




in the t. -direction. 
<P 


m , : 

<}> 


Non-dimensional applied moment at the edge. 


v = -PA 2 -. 

Y NT 2 * 

o o 


Non-dimensional surface density of the shell 



material . 
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1. INTRODUCTION 



The first paper on the dynamic plastic behavior of a 
beam was published in the early 50' s by E. H. Lee and 
P. S. Symonds [1]. In this problem and in those by sub- 
sequent authors [2/ etc.], the material is assumed to be 
rigid perfectly plastic and the vibrations of the struc- 
tures are neglected so that the problem of unloading has 
not to be considered. These assumptions, which simplify 
problems considerably, do not introduce important errors 
in many cases, since on the one hand dynamic energies are 
sometimes much larger than the elastic energy, and on the 
other hand the active duration of a pressure pulse is often 
much smaller than the corresponding natural time period of 
the structure. 

Theoretical investigations into the dynamic behavior 
of axisyrametric plates and shells were initiated with a 
paper published on circular plates by H. G. Hopkins and W. 
Prager [3]. Various boundary, geometrical and loading con- 
ditions have been considered: H. G. Hopkins and W. Prager 

solved the problem of simply supported circular plates 
under uniform dynamic load of rectangular profile [3]; A. 

J. Wang gave a solution of simply supported circular plates 
under impulsive loading [4]. Recently, M. F. Conroy pub- 
lished a paper on simply supported circular plates due to 
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dynamic circular loading [16]. The problem of clamped 

% 

circular plates was treated by A. J. Wang and H. G. Hopkins 
for impulsive load [5] and by A. L. Florence for a uni- 
form rectangular pulse load type [6]. A. L. Florence also 
solved the problem of clamped plates subjected to a 
rectangular pulse load type uniformly distributed over a 
central circular area [7]. The problem of an annular plate 
clamped along its inner edge and free along its outer one 
and subjected to impulsive load was examined by G. S. 

Shapiro [8] . 

In all these problems, the maximum displacement was 
supposed to be infinitesimal so that the strain rate- 
velocity relations are linear and only the bending moments 
are considered to be important. N. Jones examined the 
problem of finite deflection of simply supported circular 
plates under impulsive load and developed a theoretical 
procedure which retains both membrane forces and bending 
moments [9]. The solutions compare favorably with some 
experimental results on simply supported circular plates. 

Although the behavior of circular plates has been con- 
sidered in some detail, there are few solutions available 
for the axisymmetric behavior of shells. P. G. Hodge 
studied the rigid plastic behavior of an infinitely long 
cylindrical shell reinforced with equally spaced rigid rings 
uniformly loaded with a rectangular pressure profile [10]. 



17 . 



Another paper by P . G. Hodge with a numerical solution, 
examined the influence of blast characteristics on the 
final deformation of circular cylindrical shells [11] . 

In plasticity, the choice of a yield surface is of 
crucial importance. In most of the problems dealing with 
plates, the Tresca yield surface, because of its simplicity, 
has been adopted. However, for shells the Tresca yield 
surface becomes- quite complicated. Although Tresca yield 
surface for axisymmetric shells has been developed by E. T. 
Onat and W. Prager [12] , it has been used rather as a 
reference surface. To make shell problems tractable anal- 
ytically, some simplified yield surface must be devised. 

In static plasticity, the general theorems of limit analy- 
sis allow one to bracket the exact value of the collapse 
pressure and to evaluate the degree of accuracy of an 
approximate yield surface. P. G. Hodge and B. Paul have 
considered the influence of approximate yield conditions in 
the dynamic behavior of a simply supported cylindrical shell 
subjected to uniform pressure pulses having various pro- 
files [13]. Although they made some useful suggestions, no 
general theorems corresponding to the limit theorems have 
been established with respect to the predictions of approx- 
imate yield surfaces. 

The spherical shell which has two equal curvatures is 
the next shell geometry that w T as considered. R. San- 
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karanarayanan obtained the solution to the problem of 
spherical cap subjected to uniform dynamic load which decays 
exponentially [14] or has a rectangular profile [15]. The 
results obtained are involved and a numerical approach had 
to be used to examine the question of admissibility. 

It has been observed in elasticity that problems are 
usually simpler when analyzed as shallow shells and that the 
results obtained may give some indications of the correspon- 
ding behavior of deep shells. This study starts with 
dynamic plasticity of shallow spherical shells. Analytical 
procedures have been used as far as possible and three 
kinds of yield surfaces have been used, of which two only 
give distinct solutions. Then general shells have also 
been examined. The solutions obtained are complicated, for 
in the cases of plates, except for [9], and of cylindrical 
shells, there are tv/o general stresses and one velocity com- 
ponent, for which there is one equilibrium equation; the 
other two equations can be found from the yield surface and 
the flow rule. In the case of spherical shells, there are 
five general stresses and two velocity components, for which 
there are three equilibrium equations and the other four 
equations can be found from the yield surface and the flow 
rules. It is obvious that difficulties increase with the 
number of unknowns and the order of the equations. 
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2. SHALLOW SPHERICAL SHELLS 



2 . 1 General Relatio ns 

2.1.1 Strain rate-velocity relations 

For rotationally symmetrical shallow shells [17] 

loaded symmetrically, the relations between the generalized 

2 

strain rates and velocities when it is assumed that y' <<1, 



max 



y' being the slope of the meridian of the shell, are: 



v - y 'w 



x 



( 2 . 1 . 1 . 1 ) 



k = - h w ' + 
0 x 



( 2 . 1 . 1 . 2 ) 



with 



e, = v* - y"w (2.1.1. 3) 

<P 

= -h (w ' + y " v) ' (2. 1.1. 4) 

( ) • = sLLA 

[ ’ dx 



M 

and h = (2. 1.1. 5) 

AN 

o 
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In this case, the reference length A is the base radius 

* 

R of the shell: 



A = R 



Then, from (2. 1.1. 5) page 19: 

M 

. _ o 

h “ RN~ 
o 

, _ H (2. 1.1. 6) 

h 2R 

2.1.2 Equilibrium Equations 

Using the results from [17] for shallow 
symmetrical shells loaded symmetrically , ‘ and with some mod- 
ifications to include the inertia terms, the rates of 
internal energy dissipation and external energy input are: 



d. . 
int 



rl 

o (n e ®0 + v* + ■Ve + m ,j,y xax 



d 



ext 



r 1 

[x(p - yw)w + x(p, _ v)v] dx + 
n 9 

• __ • • • 

[x(sw + n^v) - hm^(w' + y"v)] x=1 



P n and p^ are respectively the normal and meridional 
components of the applied load. 

With eg, e^, kg, k^ from (2. 1. 1. 1) - (2. 1. 1. 4) page 19 
and applying the principle of virtual velocity, we obtain 
the following equations of equilibrium: 
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4 

xs = h [xm^) ' - m 0 ] - 2hp^x + -j yh 2 ft^x 

(xn^) ' - n Q - y" (xs) + p^x = xyv 
y" (xn A ) + y'n 0 + (xs) ' + p^x = xyw 



V 



n 



where = w' + y"v is the rate of slope change or angular 
velocity. has a sense of rotation contrary to the sense 

of positive m^ and the term -j yh 2 ft^ = -iC-fi^) where i is the 
non-dimensional moment of inertia of unit area of the shell 
with respect to the reference surface, represents the rotary 
inertia . 

If the load is applied only in the normal direction as 
a pressure load, and which is our case, we have: p^ = 0 , 

p = P n , and if the rotary inertia is neglected, then: 



xs = h [ (xm , ) ' - m Q ] 

9 U 



(xn^) ' - n Q - y" (xs) 



= xyv 



y"(xn^) +y'n 0 + (xs)' +px= xyw 



(2. 1.2.1) 
(2. 1.2. 2) 
(2.1.2. 3) 



For a shell without a hole, simply supported at the 
edge x = 1, and made of isotropic material, we have the 

following boundary conditions : 

At x = 0: 

From isotropy consideration: 



n 0 n cJ> 



, m 0 



From symmetry consideration: 

s = 0 



v = 0 



w ' = 0 



or there will be 
a hinge 
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At x=l : 



For simply supported edge, not free to move: 



m, = 0 
<P 

• • 

w = 0 , v = 0 

The initial conditions are: 



or there will be 
a hinge 



x = 0 , v=v = 0 , w = w= 0 



2.1.3 Yield Surfaces 

E. T. Onat and W. Prager have developed a yield 
surface for a rigid perfectly plastic material that obeys 
Tresca's yield condition [12]. The complexity of the equa- 
tion defining the yield condition is such that considerable 
difficulty may be expected in the solution of all but the 
most trivial of problems. Therefore simpler approximate 
yield surfaces have to be devised. Three approximate yield 
surfaces will be used in this study. These are: the 

two-moment limited interaction proposed by Hodge [17], the 
uncoupled square yield surface which is a simplified version 
of the former and the uncoupled diamond yield surface 
proposed by N. Jones. 
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n 



4 ) 



m 



4 ) 





Fig. 4a 



Fig. 4b 



Two-Moment Limited-Interaction Yield Surface 




m 






-s> m. 



H 



8 



Uncoupled Square Yield Surface 



Fig. 5b 





Fig. 6a Fig. 6b 

Uncoupled Diamond Yield Surface 
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The Hodge and Jones yield surfaces allow interaction 
between force and force or moment and moment but consider 
interaction between force and moment to be of limited 
importance. These yield surfaces and the planes that define 
them and the normality requirements on the associated gen- 
eralized strain rate vectors are presented in the following 
tables . 



Table 2 . 1 . 3.1 

Two-Moment Limited-Interaction Surface [ 17 ] 



Face 


Equation 


Strain Rate Vector 
(ee, e^ , £q, £9) 


1 


n e 


= 


1 


y (1,0 ,0,0) 


2 


n, 

<t> 


= 


1 


y (0,1, 0,0) 


3 


~ n e +n (p 




1 


p(-l, 1,0,0) 


4 


~ n e 


= 


1 


y .(-1,0, 0,0) 


5 


~ n t 


= 


1 


p (0 ,-l ,0 ,0 


6 


n e~ n < p 


= 


1 


y (1,-1, 0,0) 


7 


m e 


r= 


1 


p(0, 0,1,0) 


8 


m. 


r= 


1 


p(0, 0,0,1) 


9 


-m Q +m , 
y 9 


= 


1 


y (0,0, -1,1) 


0 


" m e 


= 


1 


p(0,0,-l,0) 


a 


-m. 


= 


1 


p (0 , 0 , 0 , -1) 


6 


m 0- m 9 


r= 


1 


t — 1 
1 

V. 

1 — 1 

0 

•>* 

0 

7X 
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Table 2. 1.3. 2 

Uncoupled Square Yield Surface 



Face 


Equation 


Strain Rate Vector 

<v V V V 


1 


n 9 = 1 


y (1,0, 0,0) 


2 


n , = 1 

<f> 


y (0,1, 0,0) 


3 


-n 0 = 1 


y(-l, 0,0,0) 


4 


-n, = 1 

<f> 


y (0,-1, 0,0) 


5 


m 9 = 1 


y (0,0, 1,0) 


6 


m <f, = 1 


y (0,0, 0,1) 


7 


_m e = 1 


y (0,0, -1,0) 


8 


-m , = 1 


y (0,0, 0,-1) 




Table 2. 


1.3.3 




Uncoupled Diamond 


Yield Surface 








Face 


Equation 


Strain Rate Vector 

<v V V v 


1 


V"* = 1 


y (1,1, 0,0) 


2 


— n e +n if - 1 


y (-1, i,o,o) 


3 


- n 0-% = 1 


y (-1,-1, 0,0) 


4 


n e- n <t> = 1 


y (1,-1, 0,0) 


5 


1 1 

II 

-e- 

g 

+ 

CD 

g 


y (0,0, 1,1) 


6 


1 1 

II 

g 

+ 

CD 

g 


y (0,0 ,-1,1) 


7 


- ra 9' m <J, = 1 


y (0,0, -1,-1) 


8 


me-nu = 1 


y (0,0, 1,-1) 
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Let : 

p gT = static collapse pressure of a shallow spherical 
shell structure according to the Tresca yield 
surface . 

P sH = static collapse pressure of a shallow spherical 
shell structure according to the two-moment 
limited-interaction yield surface. 

P sS = static collapse pressure according to the uncoupled 
square yield surface. 

p' sD = static collapse pressure according to the uncoupled 
diamond yield surface. 

Then according to the theorems of limit analysis 
(e.g., P. G. Hodge [18]) , we can construct the following 
bounds : 



0.618 


P sH 


1 P sT 


< 


P sH 


(2 .1. 3. 1) 


0 .309 


P sS 


± P sT 


< 


P sS 


(2.1. 3. 2) 


0.618 


P sD 


1 P sT 


< 


2p sD 


(2.1. 3.3) 



These bounds correspond to complete inscribing or cir- 
cumscribing of the yield surfaces. The interval between 
the lower and upper bounds can be reduced if only part of 
the approximate yield surface has to inscribe or circum- 
scribe the exact Tresca yield surface. It will be seen, for 
instance, that for the uncoupled diamond yield surface, the 
direct stress resultants of a shallow sphere under the 
static collapse load have a value of -1/2 approximately. 
Then, with n^ = -1/2, n^ = -1/2, the uncoupled diamond 
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yield surface can be made to circumscribe the correspond- 
ing part of the Tresca yield surface by a multiplying factor 
of 1.5 instead of 2 as started in (2. 1.3. 3) page 26. 



2.1.4 Discontinuities 

From dynamic consideration, it can be shown that 

s, n,, m, must be continuous whereas n Q , m Q may be dis- 
<j> cp y y 

continuous [17]. 



For velocities , we observe that the basic 
assumption that straight normals to the reference' surface of 
a shell remain straight and normal' to the reference surface 
during deformation implies that shear strains are negli- 
gible, and therefore w cannot be discontinuous. However 
w' and v may be discontinuous [17]. When such discontin- 
uities occur, we can see, from the expressions (2. 1.1.1)- 
(2. 1.1. 4) page 19, that the strain rate vector has the fol- 
lowing direction: 



r • % 

e 0 




r \ 

0 


• 

e , 

J * 


i. = -L • 


Av 

> 


CD 


Ax 


0 


• 




, • . 


k 4> 

V T J 




-h (Aw 1 +y " Av) 

. 



(2. 1.4.1) 



\ 

Ax 

/ 

line of discontinuity 



Ax is the increase in x to get from one side of the dis- 
continuity line to the other and Af is the jump in F cor- 
responding to Ax. 
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2 . 2 Dynamic Plastic Response of Simply Supported Shallow 
Spherical Shells made of Rigid Plastic Material that 
Obeys the Two-Moment Limited-Interaction Yield Surface 



We will solve the system of equations (2. 1.2.1)- 
(2. 1.2. 3) page 21 using successively the three yield sur- 
faces defined in (2.1.3), page 22 In each case, we begin 
by determining the static collapse pressure, both for simply 
supported and clamped supports . A solution is then devel- 
oped for the dynamic response of a shallow shell with simple 
supports when the pressure is applied dynamically with a 
peak value larger than the static collapse pressure. 

2.2.1 Static Collapse Pressures 



so that the inertia forces are negligible, the equations of 
equilibrium (2. 1.2. 2), (2. 1.2. 3) page 21 become: 



The applied dynamic pressure is assumed to have the form 



of an exponentially decaying pulse: 




or, with non-dimensional quantities: 



-t 




When the pressure is increased slowly from zero 



xs = h [fom ) * — m q ] 



( 2 . 2 . 1 . 1 ) 



(xn^) * “ n 0 - y"xs = 0 

y " (xn^) + y ‘n 0 + (xs) * + px = 0 



( 2 . 2 . 1 . 2 ) 



(2. 2. 1.3) 



The boundary conditions for a simply supported 



edge are: 




m 



<P 



x = 1 : 



0 



(2. 2. 1.4) 
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Because of the boundary condition at x=0 , a suitable 

* 

yield regime will be made of one of the faces 1, 2, 4, 5 and 
one of the faces 7, 8, 0, a. Because of the boundary con- 
dition at x=l, only faces 7 or 0 can be retained. Moreover, 
since the pressure is applied externally, we expect n^ to 
be negative and m^ to be positive, according to the sign 
convention in Fig. 3b page 19 Thus the number of choices 
reduces to two:- regime 4-7 or regime 5-7. It turns out 
that 5-7 is the correct yield regime. (Fig. 4a, b) 



With regime 5-7, we have: 
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(2. 2. 1.6) 
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With n^ = -1, (2. 2. 1.2) and (2. 2. 1.3) page 21 become: 

[ - 1 - n Q - y "xs = 0 (2. 2. 1.7) 

[y"x + y' n g + (xs) 1 + px = 0 (2. 2. 1.8) 



Solving this system and using the shallow shell approx- 
imation: y' 2 << 1, we get: 

xs = - -j px 2 + xy ' (2.2.1.10) 

n Q = - 1 + y" ( j px 2 - xy ' ) (2.2.1.11) 

By integrating (2. 1.2.1) page 21 between 0 and x, we 
obtain: 
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( 2 . 2 . 1 . 12 ) 
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A . Simply Supported Shells 
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At x=l, 




0 for simply supported edge. Thus: 
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(2.2.1.13) 



For a flat plate, y' =0, and we find again the col- 
lapse pressure of circular flat plate. 

With this expression of p, we have thus: 



0 




Fig. 7 
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- (3h + 3 



£y ' d£) x 2 + xy' 
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-1 + y " [ ( 3h + 3 



Zy 'dO x 2 - xy ' ] 



i - U + £ 



£y'd?)x 2 + ~ j Cy'dC 



(0<x<l) (2.2.1.14) 

(0<x<l) (2.2.1.15) 

( 0<x<l) (2.2.1.16) 
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The flow rules are given by (2. 2. 1.5) and (2. 2. 1.6) 
page 29. With e and k from (2. 1.1.1) and (2. 1.1. 4) 

6 (p 

page 19, we have: 

v = y’w (2.2.1.17) 

w' + y"v = A = constant (2.2.1.18) 

Solving, using the shallow shell approximation and with 
the boundary condition at x = 1, we have: 

w = w q (1 - x) (2.2.1.20) 

and v = y'w (1 - x) (2.2.1.21) 

From these velocity expressions, we have: 

e = -y 'w £ 0 since y' and w q are positive 

£ = _ ]l [_ W + yly" w (]_ _ x ) ] 

x o JJ o 

For a shallow shell, y' is small and y" which is approx- 
imately the curvature of the shell is also small, thus 
y'y"(l - x) can be considered as small compared to 1 and 
therefore, we have approximately: 

• h • • 

k = — w >0 since w >0 
x o — o 

The solution is consequently kinematially admissible. 

We assume the reference surface is smooth so that y’ exists 
everywhere, particularly at the apex x = 0. Because of 
symmetry, we have then: 
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For most practical geometries of cap, the equation of 
the middle surface can be written in the form [17]: 



For y' to exist at x = 0, we assume n > 1. 

To determine c, we use the fact that at x = 1, we must 
have y = Z/R, where Z is the total depth of the shell. 

Then : 



and the equation of the shallow shell becomes: 




n > 1 



( 2 . 2 . 1 . 22 ) 



with (y ' 2 ) = n 2 |y << 1 



max 



From (2.2.1.22), we have: 




and 



x 



£y 'd£ 




With these expressions, (2.2.1.14), (2.2.1.15), and 



(2.2.1.16) become: 



